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INTRODUCTION 


The  fail-safe  design  philosophy  requires  that 

(1)  the  damage  in  a structure  can  be  detected  before 
it  reaches  a critical  size,  and 

(2)  even  if  a critical  flaw  becomes  unstable,  the  structure 
will  have  a mode  of  crack  arrest  that  will  provide 
enough  residual  strength  and  stiffness  to  allow  time 
for  corrective  measures. 

Because  of  the  quasi-brittle  mode  of  failure  exhibited  by  graphite 
fiber  composites,  the  development  of  crack  arrest  concepts  should  be 
an  important  part  in  any  overall  structures  program  involving  these 
materials. 

The  underlying  principle  of  crack  arrestment  is  to  reduce  the 
crack  extension  (driving)  force  below  the  resisting  force  that  must 
be  exceeded  in  order  to  extend  the  crack.  This  concept  ieads  to  a 
set  of  general  classifications  of  crack  arresters,  namely: 

1.  Arresters  that  decrease  the  crack  driving  force  on  a 
propagating  crack.  (For  example,  a skin-stringer  system 
in  which  the  stringer  reduces  the  driving  force  on  the 
crack  in  the  skin.) 

2.  Arresters  that  increase  the  material  fracture  toughness 
encountered  by  the  propagating  crack.  (For  example, 
integral  "toughening"  or  "softening"  strips  in  a primary 
panel . ) 

3.  Arresters  that  both  decrease  the  crack  driving  force  and 
increase  the  material  toughness. 


CRACK  ARREST  IN  COMPOSITE  STRUCTURES 


There  are  two  distinct  basic  crack  arrest  concepts  as  described 
in  the  Introduction.  These  two  techniques  are  illustrated  in  Fig.  1. 

A primary  laminate  without  arrest  capability  is  shown  in  Fig.  la.  In 
Fig.  lb  the  crack  is  arrested  by  the  second  category  by  inserting  an 
arrester  laminate  having  a high  dynamic  fracture  toughness.  The  first 
category  is  represented  in  Fig.  lc  where  the  addition  of  a stiffener 
reduces  the  local  crack  driving  force.  In  both  Figs,  lb  and  lc  there 
is  some  return  of  the  kinetic  energy  which  must  be  “absorbed"  by  the 
arrester . 

There  are  only  a few  papers  in  the  open  literature  relating  to 
crack  arrestment  in  advanced  fiber  composite  structures.  The  buffer 
strip  (arrester  laminate;  Fig.  lb)  concept  has  been  shown  to  be  capable 
of  providing  crack  arrestment  in  graphite  fiber  epoxy  laminates. 
Eisenmann  and  Kaminski  [1]  , Bhatia  and  Verette  [2],  Sendeckyj  [3], 
and  Hess  et  al.  [4]  used  the  two-parameter  description  of  fracture  in 
composites  proposed  by  Waddoups  et  al . [5]  to  predict  the  strength  of 
panels  containing  an  arrested  crack. 

All  of  the  above  crack  arrest  analyses  for  composites  are  based 
on  static  stress  analysis.  However,  crack  arrest,  being  a dynamic 
process,  can  be  accurately  modeled  only  by  a dynamic  fracture  analysis. 
It  has  been  shown  [6]  rather  dramatically  that  statically-based  crack 
arrest  analyses  can  seriously  overestimate  a structure's  crack  arrest 
capability.  This  potentially  dangerous  statically-based  overestimation 

Numbers  in  square  brackets  indicate  references  at  the  end  of  each 
section  of  the  report. 


is  due  to  the  facts  that  static  calculations  (1)  cannot  account  for 
the  return  of  kinetic  energy  to  the  crack  tip  and  (2)  must  assume  that 
the  fracture  energy  is  a constant  that  is  independent  of  the  crack  pro- 
pagation speeds.  The  first  of  these  dynamic  requirements  is  a function 
of  the  geometry  of  the  structure  and  is  automatically  included  in  any 
fully  dynamic  fracture  analysis.  The  second  requirement  is  accommodated 
via  a dynamic  fracture  toughness  which  is  a material  property. 


DYNAMIC  FRACTURE  AND  CRACK  ARREST  ANALYSES 

If  a crack  is  extending,  the  conservation  of  mechanical  energy 
requires  that 


dW  _ dU  dT 
dA  dA  dA 


0) 


where  W is  the  work  done  on  the  structure  by  the  external  forces,  U is 
the  stored  elastic  strain  energy,  T is  the  kinetic  energy,  A is  the 
crack  area,  and  R is  the  dynamic  fracture  toughness  of  the  material. 
The  dynamic  crack  extension  force  is  defined  by 


dW^  dU  dT 
y " dA  ’ dA  ' dA  * 


(2) 


Thus,  from  the  principle  of  conservation  of  mechanical  energy, 
crack  extension  or  continued  propagation  will  occur  when 

$f(a,t)  = R (a)  (3) 

and,  no  crack  growth  occurs  or  crack  arrest  results  when 

g?(a,t)  < R (a)  (h) 

where  a is  the  crack  length,  a is  the  crack-tip  velocity  and  t is  time." 
There  is  both  theoretical  and  experimental  evidence  [7-12]  that  in  homo- 
geneous isotropic  mater ial s 1 5^  i s a function  of  crack  length  (a)  and  time 
(t) , and  that  R is  a material  property  that  is  dependent  primarily  on 

“Note  that  the  condition  ^f(a,t)>R(a)  would  violate  the  principle  of 
conservation  of  mechanical  energy  and  therefore,  cannot  exist. 
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the  propagation  velocity  (a);  hence  the  functional  dependencies  indicated 
in  eqns.  (3)  and  (A).  Whereas  R(a)  has  been  determined  for  a number  of 
homogeneous  isotropic  materials,  it  has  not  been  determined  for  any  of 
the  fiber  composites. 
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Fig.  1 Examples  of  Principal  Concepts  for  Promoting  Crack  Arrest. 

(a)  Primary  laminate  cannot  arrest  an  unstable  crack; 

(b)  Primary  laminate  with  arrester  laminate  which  increases 
the  fracture  toughness  encountered  by  crack;  and  (c)  Primary 
laminate  with  stiffener  which  decreases  the  crack  driving 
force.  (Also,  see  [6].) 
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SECTION  2 

APPROXIMATE  MODE  I VELOCITY  CORRECTION  FACTOR  FOR 
0°  UNIDIRECTIONAL  FIBER  COMPOSITES 


ABSTRACT 

An  approximate  Mode  I dynamic  crack  extension  force  for  a 
highly  orthotropic  0°  unidirectional  fiber  composite  is  derived  in 
terms  of  the  quasi-static  crack  extension  force  and  a velocity  correc- 
tion factor.  The  analysis  is  based  on  the  dynamic  extension  of  a 
static  Mode  III  i sot ropi c/Mode  I orthotropic  analogy  which  was  first 
recognized  by  F.A.  McClintock.  The  complete  dynamic  analogies  are 
given  and  the  derived  orthotropic  velocity  correction  factor  is  com- 
pared with  various  isotropic  velocity  correction  factors. 
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INTRODUCTION 


> 


The  use  of  linear  elastic  fracture  mechanics  (LEFM)  in  the  study 

of  the  onset  of  crack  instability  in  brittle  homogeneous  isotropic  materials 

is  one  of  the  modern  accomplishments  in  mechanics.  The  onset  of  rapid 

fracture  may  be  characterized  by  critical  values  of  the  static  stress 

intensity  factor  or  the  quasi-static  crack  extension  force  gf  which 
have  been  shown  [I]  to  be  related,  for  example,  as 


for  Mode  I I I 


where  G is  the  shear  modulus  of  the  material.  (Similar  expressions  can 
be  written  for  Modes  I and  II.) 

At  present,  it  is  generally  accepted  that  for  brittle  homogeneous 
isotropic  materials,  the  dynamic  crack  extension  force ^ on  a crack 
moving  at  constant  velocity  can  be  expressed  as  the  product  of  ^ and  a 
velocity  correction  factor  g(v)  as  [2,3] 


S5  ■ 9<v)g? 


where  g(v)  depends  on  the  fracture  mode  and  is  a function  of  the  material's 
di latational , shear  and  Rayleigh  wave  speeds. 


rePresents  the  nonrecoverable  loss  of  energy  per  unit  of  crack 

extension  and  is  also  called  the  strain  energy  release  rate.  We  prefer 
to  call  it  the  crack  extension  force  or  the  crack  driving  force  which 
we  feel  is  particularly  appropriate  for  dynamic  analyses. 


The  purpose  of  this  section  is  to  present  the  derivation  of  an 


r 


i 

i 


approximate  velocity  correction  factor  w(v)  which  is  applicable  to  Mode  1 
dynamic  crack  propagation  in  0°  un id i rect Eona I graphite  fiber  composites. 
The  fiber  composite  is  modeled  as  an  equivalent  homogeneous  orthotropic 
material.  The  analysis  is  based  on  an  extension  of  an  analogy  between 
the  static  Mode  III  crack-tip  stress  and  displacement  fields  in  homogeneous 
isotropic  materials  and  the  static  Mode  I crack-tip  stress  and  displacement 
fields  in  homogeneous  highly  orthotropic  materials.  This  analogy  was  first 
observed  by  McClintock  [4].  Therefore,  to  the  extent  that  it  is  useful  to 
define  a crack  extension  force  or  a critical  stress  intensity  factor  for 
fiber  composite  materials  based  upon  such  stress  and  displacement  fields, 
the  dynamic  crack  propagation  in  a fiber  composite  may  be  described  in 
general  terms  by 

S£--(v)^  (3) 

c c 

where  and  & are  the  dynamic  and  static  composite  crack  extension 

forces,  respectively. 


THEORETICAL  ANALYSIS 


Crack-Tip  Fields 

In  an  unpublished  paper,  McClintock  [ 4 ] described  an  analogy 
between  the  Mode  III  isotropic  and  the  Mode  I highly  orthotropic  crack- 
tip  stress  and  displacement  fields.  Both  the  Mode  III  isotropic  case  and 
the  Mode  I highly  orthotropic  case  are  shown  in  Fig.  1.  The  Mode  I highly 
orthotropic  equations  are  approximate  and  are  based  on  the  assumption  that 
E2  > > E j , where  E^  and  Ej  are  the  extensional  moduli  orthogonal  to  and 
along  the  crack  direction,  respectively.  This  highly  orthotropic  require- 
ment is  generally  met  by  0°  unidirectional  graphite  fiber  epoxy  composites. 

The  static  analogy  provides  functional  equivalence  between  the 
following  variables: 


Mode  III  I sotropic 


Mode  I Orthotropic 


■t  *r,  * t 


where  the  primed  variables  refer  to  Mode  III  isotropic  and  the  unprimed 


variables  refer  to  Mode  I orthotropic.  The  x.  are  the  coordinates  with 
the  associated  u.  displacements  and  the  a.j  stresses.  G is  the  isotropic 
shear  modulus  and  G ^ 2 's  the  in-plane  orthotropic  shear  modulus.  The 
static  stress  intensity  factors  are  the  k.  which  are  defined  in  terms 
of  the  stresses  at  infinity  a..  and  the  crack  lengths  c1  and  c.  A 
complete  derivation  of  this  analogy  is  given  in  the  Appendix. 

The  extents  to  which  the  relevant  equilibrium  and  constitutive 
equations  are  satisfied  by  the  analogy  will  be  examined  next.  These  are 


Mode  III  Isotropic 


3°i  1 a 1 3o» 1 a 1 

Equilibrium:  3 + 3 — = 0 

dx  j , <*x2 , 


or. 


Mode  I Orthotropic  (Via  Analogy) 

_ a(A7^r°i2)i  *°22 

4/G|2/E2  *2) 


3x, 


= 0 


3c  3^22 

— — + -r-==-  “ 0 
3Xj  3x2 


(5) 


Stress-Strain:  a.,,,  = G 


3u3, 


1 *3*  3x,, 


3u_ 

°12>/E2/G12  * >/E2G12  STj" 


or, 


3u2 

°12  " G12  3x. 


(6) 


3u3, 


3u„ 


Stress-Strain:  aol-,  * G r-  - 

J , 


or, 


3uj 

°22  " E2  3x1 


(7) 
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The  equilibrium  and  constitutive  equations  for  a 0°  unidirectional  com- 
posite as  shown  in  Fig.  1(b)  are  [5] 


Equi I ibrium: 


Stress-Strain: 


> -o 

3xj  3x2 


, . „ + 

12  ^I2^3x2  Sx,/ 


°22  " V2I°11  " E2  3x2  * 


It  is  clear  that  eqns.  (5)  and  (8)  are  equivalent.  However,  the  equi- 
valences between  eqns.  (6)  and  (9)  and  between  eqns.  (7)  and  (10)  are  not 

so  apparent.  Eqns.  (6)  and  (9)  would  be  approximately  equivalent  if 
3u.  3u2 

•5 — may  be  neglected  when  compared  with  -5 — • In  the  Appendix  this 
ox 2 

is  shown  to  be  consistent  with  other  approximations  which  have  already 
been  made  in  deriving  the  analogy,  eqns.  (4).  Also,  eqns.  (7)  and  (10) 
would  be  approximately  equivalent  if  a22  » v2lall  ‘ This  condition  is 
shown  to  be  satisfied  also  in  the  Appendix.  Thus,  the  governing  orthotropic 
field  equations  are  consistent  with  the  respective  equations  obtained  via 
the  analogy. 

Extension  of  Analogy  into  the  Dynamic  Regime 

In  order  to  investigate  the  dynamic  fracture  of  composites  via 
the  analogy,  it  is  necessary  to  extend  it  into  the  dynamic  regime.  The 
additional  analogous  variables  which  must  be  obta ined  a re  the  time  and  the 
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crack-tip  velocity.  (To  date  the  crack-tip  acceleration  has  not  entered 
the  solutions  for  the  dynamic  crack  driving  forces  in  isotropic  materials, 
and  thus  the  analogies  developed  here  will  not  consider  acceleration. 
However,  we  reserve  judgment  as  to  whether  or  not  acceleration  should  be 
included  in  dynamic  fracture  analyses,  and  we  shall  consider  this  issue 
in  a subsequent  report.) 


1 . Analogy  Between  Times 


The  equations  of 

motion 

for  Mode 

III  i 

isotropic  and  Mode 

1 ortho- 

tropic  bodies  are 

So,,,, 

3o„ 

32u,, 

Mode  1 1 1 Isotropic: 

2*3 

— r\  1 

, 3 ' 

(11) 

3X|  , 

3x2, 

P 

3t 1 2 

Mode  1 Orthotropic: 

3*.2 

3x, 

3<t22 

3x2 

= P 

32u2 

3t2 

(12) 

where  p'  and  p are  the  Mode  III 

isotropic 

and 

Mode  1 orthotropic 

mass 

densities,  respectively, 

and  t ' 

and  t are 

the 

respective  times. 

Assumi ng 

that  there  is  an  analogy  between  t'  and  t such  that  t'  -*■  at , where  a is 
a proportionality  constant,  a may  be  determined  by  substituting  the 
analogy  represented  by  eqns.  (1*)  into  eqn.  (11)  which  gives 


3a 


22 


s(^7T2  x2) 


sV 


3x, 


3 (at) : 


or. 


3a 


12 


3a 


22 


3x, 


?x„ 


VE2  a2  3t2 


(13) 


J 


i 
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In  order  that  eqns.  (12)  and  (13)  are  equal,  it  is  required  that 


Therefore,  the  analogy  between  times  is 


Mode  111  I sotropic 


Mode  I Orthotropic 


GW'- 


2 . Analogy  Between  Velocities 


Assuming  a constant  crack-tip  velocity,  the  Mode  III  isotropic 


velocity  V1  may  be  found  from 


V'-f  = X,,  - 


where  Xj  ■ 's  the  horizontal  axis  (xj  , direction)  in  a fixed  reference 
frame  and  c^  is  the  initial  (or  reference)  crack  length.  Also,  assuming 
that  there  is  an  analogy  between  the  Mode  III  isotropic  and  the  Mode  I 
orthotropic  velocities  such  that  V'  -*•  0V,  where  (3  is  a proportionality 
constant  and  V is  the  Mode  I orthotropic  velocity,  eqns.  W , (15)  and 


(16)  give 


-X|-cc 
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where  ^ is  the  Mode  I orthotropic  horizontal  axis  in  a fixed  reference 
frame  and  cq  is  the  initial  (or  reference)  crack  length.  Furthermore, 
from  kinematics 


V*t 


X, 


-c 


(18) 


Combining  eqns.  (17)  and  (l 8)  gives 


Therefore,  the  analogy  between  the  crack  propagation  velocities  is 


Mode  III  Isotropic 


Mode  I Orthotropic 


V' 


(20) 


3.  Analogy  Between  Crack  Driving  Forces 


Kostrov  and  Nikitin  [6]  proved  that  for  an  elastic  body  in  the 
>ence  of  any  thermal  effects  on  the  fracture  process 


0)  = -v  1 
e 


. f 

im  I 

-*-o  Jx 


1A 

IB 


' uj,i) 


x2  = +e 


dx. 


(21) 


X2  = -E 


^Note  that  by  considering  the  differential 
this  derivation  could  have  been  performed 


forms  of  eqns.  (16)  and  ( 1 8 ) , 
for  arbitrary  crack  velocities. 
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where  w is  the  power  of  energy  of  fracture  per  unit  length  of  crack, 

v is  the  crack-tip  velocity,  and  ( ) . denotes  differentiation  with 

» ' 

respect  to  Xj  . The  integration  contour  for  eqn.  (21)  is  shown  in  Fig.  2. 
Eqn.  (12)  is  valid  for  arbitrarily  changing  velocity  as  well  as  for  con- 
stant velocity. 

For  the  Mode  III  isotropic  body,  the  only  nonvanishing  stress 
component  that  contributes  to  w is  the  shear  stress  • Then,  for 

this  case,  eqn.  (21)  can  be  written  as 


“l  I I (iso) 


-V 


x2,  = +e' 


(22) 


where  aJ|||(;so)  's  Mode  III  isotropic  power  of  energy  of  fracture 

per  unit  length  of  crack.  Also,  for  the  Mode  I highly  orthotropic  body, 

the  only  nonvanishing  stress  component  that  contributes  to  u)  is  the 

i* 

normal  stress  a ^ . So, 


(23) 


where  u)j  (ortho)  ’S  t^1e  Moc*e  * orthotropic  power  of  energy  of  fracture 
per  unit  length  of  crack. 


^Due  to  Mastic  symmetry  with  respect  to  the  Xj  axis,  the  shear  stress 
Oj2  along  x^  is  zero  and  does  not  contribute  to  w.  See  eqn.  (A22)  in 
the  Appendix. 
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The  relationship  between  uk  and  the  dynamic  crack  extension 

force  <&..  is  given  by  [6] 
a . j 


“j  " 


W) 


So,  from  eqns.  (22)  and  (2*») , the  Mode  III  isotropic  dynamic  crack 
extension  force  ^;||)(Jso)  is 


d : 1 1 1 ( i so) 


""o  Lv,  (0^'3'^v) 


xi  ■ +e' 


dx, , 


(25) 


X * * -£ : 
x2  E 


From  eqns.  (23)  and  (2**) , the  Mode  I orthotropic  dynamic  crack  extension 
force  K.  I (ortho)  is 


1 (ortho) 


-1  im 
e-*-o 


a,b  \ i / 


x„  ■ +e 


dx. 


(26) 


x-  = -e 


It  is  now  necessary  to  establish  the  analogy  between 

and  I (ortho)  * By  usin9  the  analo9y  for  a2'3'  * u3 ' * x 
which  is  given  in  eqns.  (4) , eqn.  (25)  gives 


1 1 1 (isc) 
j,  and  x2, 


-2k- 


I 


I 


I 


i 

! 


where  ^.^naj  represents  the  Mode  III  i sot ropi c/Mode  I orthotropic 
analogous  dynamic  crack  extension  force.  Egns.  (26)  and  (27)  differ 


FvT 

only  in  the  small  quantities  £ and  — • e in  the  respective  equations. 
However,  because  e is  arbitrary  and  tends  to  zero  in  the  limit,  eqns.  (26) 
and  (27)  are  equivalent  and  therefore  ^j.|(0rt^0)  's  analogous  to 
^d- 1 II  (iso)  * Furthermore,  because  of  the  direct  analogy^  between 

the  as  represented  by  eqns.  (25)  and  (26),  a direct  analogy 

between  the  & . . follows  because  in  the  limit  as  v -+•  0,  eqns.  (25) 
s : j 

and  (26)  represent  the  expressions  for  the  & . Thus, 

s : j 


Mode  111  Isotropic 


Mode  I Orthotropic 


s: 1 1 I (iso) 


1 1 I (iso) 


^s : I (or 


(ortho) 


d: I (ortho) 


(28) 


^The  phrase  "direct  analogy"  is  used  to  indicated  that  no  multiplicative 
constants  are  required  in  order  to  complete  the  analogy.  For  example, 
the  analogy  between  ar|d  *2  's  not  "direct"  since  the  multiplicative 

constant  ^/G ^ 2 / E 2 is  required  to  complete  it. 
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f < 


VELOCITY  CORRECTION  FACTOR 


By  using  the  analogies  which  have  been  presented,  a velocity 
correction  factor  may  be  obtained.  The  velocity  correction  factor  for 
the  Mode  III  isotropic  case  is  [8] 


9m(V> 


. V' 

I + — 
c 

s 


1/ 

where  cg  = (G/p1)  2 * * * * * is  the  shear  wave  speed  and  p‘  is  the  mass  density 
of  the  isotropic  medium. 

By  assuming  eqn.  (2)  to  be  a concept  which  is  generic  to  dynamic 
fracture  in  any  material,  the  following  equations  may  be  written: 


^d:  I II  (iso)  9 1 1 1 ^ ^s:lll(iso) 


I (ortho)  " w|  (0°) (V)  ^s:  I (ortho) 


I 

I 

I 


where  the  Mode  I orthotropic  velocity  correction  factor  w| (q°)  has  keen 

written  with  subscripts  to  denote  its  relationship  to  Mode  I fracture  in 

0°  unidirectional  composites.  Because  of  the  direct  analogies  between  the 

and  between  the  as  stated  in  eqns.  (28),  g(||  and  w( 

are  directly  analogous  also.  Thus,  utilizing  eqns.  (4),  (20)  and  (29) 

provides 
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WI(0°)(V) 


(E2/G12)V>  V 
c ) '' 


(E2/Gi2  r V 


I2E2/Pc  ) 


w| (0°) 


/ \V 

where  Cs  = ^ Pc ) 2 's  the  composite  shear  wave  speed  and  pc  is  the 

composite  mass  density. 

Eqn.  (31)  may  be  expressed  in  terms  of  the  longitudinal  wave 
speeds  also.  For  example,  along  the  stiffer  modulus  direction  (x»)  the 
longitudinal  wave  speed  is  given  by  C 2 = , and  using  this 

definition  in  Cg  gives 


(if)' 


Substitution  of  eqn.  (32)  into  eqn.  (31)  gives 


wi(o°)(v)  = — 


-Mi 

< • ft)*  4 


. * 1 — •• 


---.v  IM 


Eqn.  (31)  is  plotted  in  Fig.  3.  For  comparison,  the  Mode  I 


k 


velocity  correction  factors  for  homogeneous  isotropic  elastic  materials 

as  calculated  by  several  researchers  are  shown  also.  The  Broberg  (9) 

2 

result  which  is  for  Poisson's  ratio  v = y is  for  a crack  which  at  time 
t = 0 grows  symmetrically  from  zero  length  at  constant  velocity.  The 
Freund  [10]  result  which  is  for  V = is  for  a semi - i nf in i te  crack 
which  at  time  t = 0 grows  at  constant  velocity.  The  Rose  [8]  result  is 
simply  a linear  approximation  of  previous  results.  The  most 
striking  feature  of  the  w| (q°)  computed  here  is  that  it  does  not  go 

V 

to  zero  at  ~0.92  which  represents  the  Rayleigh  wave  speed  for  homo- 
s 

geneous  isotropic  materials.  Also  for  V greater  than  zero,  w| (0°)  is 
always  greater  than  the  isotropic  velocity  correction  factors. 


1 


* 


SUMMARY 


The  complete  results  for  the  Mode  III  i sotropic/Mode  I ortho- 
tropic dynamic  crack  propagation  analogy  are  summarized  in  the  Table 
below. 


TABLE:  Summary  of  Mode  III  Isotropic/Mode  I Orthotropic  Dynamic 

Fracture  Analogy 


Mode  III  Isotropic 


X1 1 * (c ' ) 


Mode  I Orthotropic 


Xj  , (c) 

x2>/G12/E2 


2 ' 3 ' 

k3  ' = a2 ' 3 '« 


fs: 1 1 1 (iso) 


'd:  1 1 I (iso) 


g,M(V) 


o,2  n/E2/G1 2 

a 

22 

kl  " ° 22°°  ^ 

(G, 2/E2 ) ^ t 
v 

% : I (ortho) 
85  : I (ortho) 
W|(0°)(V) 


This  analogy  is  valid  under  the  assumptions  that  homogeneity  and  high 
orthotropy  are  adequately  represented  in  the  composite. 
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(a) 


a)  Mode  III  isotropic  case. 

b)  Mode  I 0°  orthotropic  case. 


Mode  I 


O.l  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 

V 

Normalized  Crack-tip  Velocity,  -pr- 


Fig.  3 Mode  I velocity  correction  factors  for  homogeneous  isotropic 
elastic  materials  (g J ) and  for  strongly  orthotropic  0°  uni- 
directional composites  (w|(o0)^" 
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APPENDIX 


The  purpose  of  this  appendix  is  to  present  the  results  of 
McClintock  in  greater  detail  than  that  which  is  given  in  [I*] . Although 
the  essential  ideas  and  results  are  the  same  as  those  derived  by  McClintock 
and  thus  the  credit  for  this  analogy  is  his,  we  believe  that  the  detailed 
presentation  given  herein  provides  additional  delineation  of  the  appli- 
cability of  the  analogy. 

Order-of-Magn ? tude  Approximations  of  Mode  I Orthotropic  (0°)  Elastic 
Elastic  Parameters 

The  strain  (e,y)  - stress  (a)  relations  in  principal  material 
coordinates  for  an  orthotropic  lamina  under  plane  stress  are  [5] 


'll 


'22 


12 


11 


a, 2 0 


a2,  a22  0 


66 


'11 


22 


'12. 


(Al) 


where 


11 


al 2 = a21 


22 


'12 


21 


66 


’12 


Ej  and  are  the  extensional  moduli  in  the  Xj  and  x^  directions, 
respectively,  is  the  shear  modulus  in  the  Xj-x^  plane,  and  the 
V.  . are  the  Poisson's  ratios  for  transverse  strain  in  the  j direction 
due  to  stress  in  the  i direction.  The  Mode  I orthotropic  stress  and 
displacement  fields^  are  given  in  terms  of  the  roots  p.  of  the  follow- 
ing equation  [7] : 


- (2., 


2 * a66/lJ  * a; 


which  has  the  solutions 


2 2 

- 


(U\2*a6t)  i (2a!  2 + a6&)  V’ 


*alla22 

(2*12*aJ2 


* For  a 0°  graphite  fiber  composite  (fibers  extending  along  the  direction), 
typical  constitutive  parameters  give 


*2 

■==■  a 20 
E1 


V21*  J 


Some  simple  combinations  of  the  elastic  constants  will  be  eval- 
uated for  the  purpose  of  order-of-magn i tude  approximations.  Using 
eqns.  (A 1 ) and  (aA) 


rThe  explicit  equations  for  these  fields  will  be  given  later  in  this 
Appendix.  -40- 


Hence,  from  eqns.  (A5)  and  (A6) , it  follows  that 


a66  ^ 2 1 ai 2 ^ and 

(2aI2  + a66 )2  >> 

Sla22  . 

Observing  that  JT- a~  can  be 

approximated  by 

1 - y when  0 < a « 

eqn.  (A3)  can  be  written  as 

” (2a  1 2 + a66  ) - (2a12  + a66) 

2 2 

1 1 * ii*  ■ ■ 

r,  i 4ana22  ] 
- 2 (2a12*a66)\ 

**  1 > ^2 

2all 

which  gives 

y2  ~2(2a 1 2 + a66 )2  + 2alla22 

' 2al  1 (2a  1 2 + a66  ) 


(2al 2 + a66) 
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By  using  the  first  condition  in  ( A7 ) on  y2  an<^  both  conditions  in  (A 7 ) 
on  y2^  , eqns.  (A8)  can  be  reduced  further  as 


1 1 ^ ^ 

y,  ~ 


2an  + a66 

an 


i?  » - 


22 


2a12  + a66 


a66 


a22 

a66 


(A9) 


Also,  the  ratio  of  these  roots  is 


v/8cT  ~ 9 . 


(A10) 


as 


Next,  consider  the  complex  variables  Zj  and  z^  which  are  defined 


Zj  ® cosG  + y^sinG 

(All) 

z2  = cosG  + y^sinG 

where  0 is  the  counterclockwise  angle  measured  from  the  x^  axis  to  r,  thus 
forming  an  r,  0 crack-tip  polar  coordinate  system.  Because  y^  and  y2  are 
pure  imaginary  quantities  (Refer  to  eqn.  (A9).),  eqns.  (All)  are 

Zj  = cosG  + i | vi | | sinG 

(A  12) 


cosG  + i |y, | sin0 


which  can  be  written  as 


Zj  » [cos20  + I M J I 2 s i n20  ] 


% « 4> , 


r i l/i 

z2  * [cos20  + I u2 1 2 sin20  J e 


(A  13) 


where  the  <f>.  are 

ri 


= tan 


.j  |P||  sin0 


COS0 


So, 


I 

I 

I 


* 


f I 


1 

.Vi 


1 4 (lv>j  I2  - l)sin20 
i + (|U2|2  " l)sin20 


1/4 


which  has  a maximum  when  ‘t’l  2 = T anc*  ® = T’  anc*  a m'n'mum  when  <f>j  ^ = 0C 
and  0=0°.  Thus, 


(A14) 


and 


(A15) 


min 


Combining  eqns.  (AI0)  and  (A  15)  gives  the  following  maximum  and  minimum 
ratios: 
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r * . % t '^+1  •>  * »• 


i 


The  singular  crack-tip  stress  field  is  given  by  [7] 


i 


I 

I 

I 

I 


I 

I 


(A  18) 


where  kj  is  the  stress  intensity  factor  and  Re  denotes  the  real  part  of  a 
complex  function. 


Recal 1 ing  that 


is  maximum  (~9)  for  6 = 0°  and  minimum 


(~3)  for  9 = j and  by  letting  y^  - y^  w y^  in  accordance  with  eqn.  (A  1 0) 
the  first  of  eqns.  (A  18)  may  be  written  as 


Re 


(A  1 9) 


-Z}5- 


. 


m i A 


V V 

where  has  been  dropped  in  comparison  with  Pj/Zj  . Substituting 

eqns.  (A9)  and  (All)  into  eqn.  (A  1 9)  gives 


a 


n 


k*  Re 

'(-1)1 

i \ 

( M 

! >/a66/all  i 

liv/a22/a66  1 

v6T 

_ 

cos©  + iya^g/ajj  sin0 

V’J 

or 


11 


Re 


>/a22/all 

(x,  + Va66/aj  i" 


(A20) 


where  Xj  = r cos0  and  = r sin0  have  been  used.  It  is  important  to  note 

1/  i / 

that  in  deriving  eqn.  (A20)  , the  requirement  that  Vt]^z\>>^2^Z2  ^as  1366,1  usec1' 
This  is  true  as  0 + 0°.  (As  0 -*■  j , p^/z^2  wh'ch  would  result 

in  about  a 25%  error  in  eqn.  (A20)  at  0 = j .)  It  is  important  to  note 
that  (1)  Oj j is  not  used  in  the  analogy  for  the  derivation  of  the  crack 
extension  force  and  (2)  for  self-similar  Mode  I crack  extension,  the  value 
of  <J|j  at  0 = 0°  would  be  the  important  value  of  j if  j were  used 
in  a given  analysis. 

Recalling  the  second  of  eqns.  (A  1 7)  and  letting  p^  -p^ssPj  , the 
second  of  eqns.  (A 1 8)  may  be  written  as 


a 


22 


or  using  eqns.  (A9)  and  (All), 


-'(6- 


ic  (0°)  Composite 


The  crack-tip  displacement  field  is  given  by  [7] 


where 


Pj  = any.  + aJ2  - aj6y. 


ql  " a12Mi  + y.2  _ a26 


where  a^  = a^  - 0 for  a 0°  orthotropic  composite.  By  letting  y^  -y2«yj  , 
eqns.  (A23)  become 


Uj  w k(>/5r"  Re 


u2  ss  k yjl7  Re 


'pi  uT  V + p’2 


V*2  1/  1/ 

-q.  — z,  + q,z/2 
w y,  1 H2  2 


(A24) 


Examination  of  the  first  of  eqns.  (A2l*)  gives 


-p'*7 


[allyi  + a12](^)%[a66  " al2]/ 


/alia22 


*66 


where  eqns.  (A9)  have  been  used  for  the  y.  . From  the  first  condition 


(A?)  ’ a66  >> I a 1 2 1 » "pl  (yf)25/ 


in  eqn. 


alla22  ’ 


(A25) 


Further,  using  eqns.  (A9)  again  gives 


P2*al2 


ana22 

a66 


. 1 1 
60E  | 


(A26) 
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where  the  definitions  of  the  a. ^ (eqns.  (A  1 ) ) and  the  typical  constitutive 
parameters  (eqns.  (A1*))  have  been  used.  The  first  terms  (a  12  ’ ' 6^)  °n 
each  side  of  eqn.  (A26)  correspond  to  each  other  and  the  second  terms 
/_  a j | a 22  1 \ 

^ — , “ Ifoi- y on  each  side  of  eqn.  (A26)  correspond  to  each  other. 

Thus,  as  eqn.  (A26)  indicates,  a^  is  not  sufficiently  larger  than  a ^ j a22^a66 
to  allow  neglecting  a j j a22^a66  'n  eqn’  ’ Theref°re»  substitution  of 

eqns.  (A25)  and  (A26)  into  (A24)  gives 


Re 


1 1 a22  (X,  +M1 X2  ) 2 + (a12  ‘ a66  )(X1+U2X2  ) 


(A27) 


Examination  of  the  second  of  eqns.  (A24 ) gives 


V»  / 2 

V2Zl  = »a”y'  + a 


1 2M1  22 


)(?)•;*.(•»  -=e)/. 


a!2a66W!n!22  2v2  (A28) 


1 


where  eqns.  (A9)  have  been  used  for  the  y..  From  the  definitions  of  the  a.j 
(eqns.  (Al))  and  the  typical  constitutive  parameters  (eqns.  (A4))  note  that 


a22  ~ W 


and 


12 


a12a66  _ 1 
17j SOG 


12 


and  therefore. 


(A29) 


Substitution  of  eqn.  (A29)  into  eqn.  (A28)  gives 


W,  ~ 2 “12^  • 


Also, 


(A3 


v2 

VlZ2 


(A3 


where  eqns.  (A9)  have  used  and  the  condition  a,,  »|a,J  has  been  used. 

oo  1 2 1 

The  ratio  of  eqn.  (A3 1 ) to  eqn.  (A30)  gives 


‘loW.z 


V2 


2H  2 


■wr 


tr  -©  £ RK)(-“8,(t  *•') 


(A3 


where  eqns.  (A  1 ) and  (A**)  have  been  used  to  evaluate  ^a^/a^j  ar>d  eqns 
(A10)  and  (A  1 5)  have  been  used  to  evaluate  (z^/z^  . So, 


Va 

W2 


■qi¥i 


Va 


a 16  to  A8. 


(A33 


and  thus. 


W2 


Va 


» 


qi¥i 


Va 


Hence,  the  second  of  eqns.  (A24)  may  be  written  as 


term  has  been  dropped. 


1/ 

where  in  accordance  with  eqn.  (A3 3 ) the  q^^z'2 
Therefore,  by  using  eqns.  (A9) , eqn.  (A3*0  becomes 


u 


2 


Re 


where  as  before,  Xj  = r cos9  and  = r sin0  have  been  used. 


Analogy  Between  Mode  III  Isotropic  (Exact)  Fields  and  Mode  I Orthotro 
(Approximate)  Fields 


The  crack-tip  stress  and  displacement  fields  for  the  Mode  III 
isotropic  case  are  [7] 


°1  '3' 


"k3 ' 


—j==r  Re  ' 

2 L(x,,  + ix2()V2  J 


a2'3' 


_ 3 


R-j  i 


Re  

^ L(x11  + ix21)V2  -I 


Uj,  = kyj2  Re  [-i  ^ (x,  , + ix2t)  /2] 


where  the  primed  variables  refer  to  the  Mode  III  isotropic  analysis, 
paired  variables  to  be  investigated  are 

Oj2  and  Ojiji  £eqn.  (A22)  and  eqn.  (A36-a)j 


(A35) 


(a) 


(b) (A36) 


(c) 

The 
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ana  a 


j^eqn.  and  eqn.  (A3b-b)J 


u2  and  u^,  ^eqn.  (A3 5 ) and  eqn.  (A36-c) J . 
Rewriting  eqns.  (A22)  and  (A36-a)  gives 


kl 

0,2  ~yT  Re 


i y a_9/a 


22  66 


1 v/ra22/a66 


[(xl  + 5 v/a66/a,  i x2)  /2  (*,  + Va22/a6  xj  /z  j 


and 


a 


1 '3' 


+ 


(A3< 


which  indicate  that  if  the  first  set  of  terms  in  the  square  brackets  in 

eqn.  (A22)  were  not  present,  a direct  analogy  between  (A22)  and  (A36-a) 
would  result. 

As  Xj+0  for  x2*0,  the  first  term  in  eqn.  (A22)  is  smaller  than 
the  second  term,  and  their  ratio  is 


where  eqns.  (Al)  and  (A*0  have  been  used.  As  -*•  0 for  Xj  > 0,  (i.e.f  in 
front  of  the  crack  tip),  both  terms  vanish  and  the  difference  is  not  impor- 
tant. As  x.  •+■  0 for  x.  < 0,  (i.e.,  along  the  crack  flanks),  the  shear  stress 
i.  i 

vanishes  for  all  jjj  . 

As  X£  -*■  0 for  Xj  > 0,  both  and  Ojijt  tend  to  zero  in  the  same 

functional  manner.  So  qualitatively,  each  term  of  eqn.  (A22)  behaves  the 

ame 

I 1 3 1 does  along  the  horizontal  axis  for  x.  > 0.  Because  the  second 
term  is  three  times  larger  than  the  first  term  in  eqn.  (A22) , the  first  term 
will  be  neglected.  However,  as  McClintock  [4]  observed,  the  neglection  of 
the  first  term  may  not  be  justifiable,  in  general,  because  the  first  term 
reduces  the  magnitude  of  the  shear  stress  above  and  below  the  crack  tip, 
where  the  shear  stress  is  maximum,  by  a factor  of  about  two.  Despite  this 
general  reservation  regarding  the  shear  stress  a ^ , it  poses  no  difficulty 
to  the  analysis  presented  in  the  body  of  this  report  because  and  U£ 
along  the  horizontal  axis  are  the  only  analogous  quantities  required.  The 
reduced  forms  of  eqn.  (A22)  and  eqn.  (A36-a)  may  be  written  respectively,  as 


where  c'  and  c are  the  respective  Mode  III  isotropic  and  Mode  I orthotropic 


crack  lengths.  Thus,  c' — ►c  represents  a useful  analogy  which  follows 
directly  from  eqns.  (A40) . (This  direct  analogy  between  c'  and  c may  be 
drawn  from  the  direct  analogy  between  Xj , and  Xj  also  since  the  respective 
crack  lengths  simply  represent  specific  lengths  along  the  horizontal  axes.) 
And,  for  eqns.  (A35)  and  (A36-c) 


u2  « kp/2"  Re  {-i  yafifia„  x,  + 


66  22  ll 


•(m  -or 


u j 


-i  i (x,,  + i*2,)V’ 


(A35) 


(A36-c) 


from  which  the  following  analogies  can  be  drawn: 


'a66a22 


1 2^2 


(A41 ) 


It  may  be  noted  that  eqns.  (A39) , (A^O)  and  (A4l)  are  consistent  with  each 
other  including  the  definitions  for  kj  and  k^ , which  when  expanded  give 


’ *.  «*  4 V ’ft 


Further  Comments 


1.  Limit  of  Approximations  - 


The  initial  level  of  approximations  which  were  used 


and  y* 


a66  1 


(2  a66 


can  not  be  relaxed  if  the  Mode  III  i sot ropi c/Mode  I 


orthotropic  analogy  is  to  remain  applicable.  If  the  more  accurate,  but  still 

2a , j ^ 3/ / a^o 

approximate,  expressions  y*  ~ - and  yj?  « - - had 

1 all  2 2a12  + a66 

been  used  (See  eqn.  (A9).),eqns.  (A21),  (A22)  and  (A35)  would  have  become, 
respectively, 


* 4 


3u.  i 

If  a check  is  made  as  to  whether  or  not  a,,,,  = G --  ■ results 

2 3 3x2, 

I 9u2  t 

by  analogy,  the  following  occurs: 

a22  <jx2 


in 


'22 


‘66 


3u„ 


a22a66 


9 Iva22/^2al2  + a66^ 


or. 


J2a 12  * a66  l2a 12  + a66  9u2 

a22a66  a22  9x2 


°22  " 


2a12  + a66 
a66 


J_  # ^2 
a22  9x2 


3u„ 


(Al»3) 


22  _ a 3^  » • «■  is  necessary  to  invoke 


I 0 

In  order  to  obtain  the  equation  = y-=-  , it  i 

a22  9x2 

the  further  approximation  that  a&&  » 2 | a 1 2 | , which  is  what  has  been  used 
throughout  this  appendix. 


2.  Mode  I Orthotropic  a)2  - y)2  Approximation  - 


In  order  to  insure  the  satisfaction  of  the  a]2~*]2  ecluat'on  in 
accordance  with  the  analogy,  the  constitutive  equation 


/3u,  3u2\ 

a,2‘  + 


This  is  done  in  the  main  text  of  this  report. 
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must  consistently  reduce  to 
appendix  that 


0,2  Ui  /■ 


Recall  from  earlier  in  this 
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' ^a22/a66 


'^a66/aH 


iv/a22/a66 

^a22/a66  x2) 1/2  - 


and 


^ '3'~"y2  Re 


i 

.(xl 

1 + fx2')l/2  . 

(A36-a; 


and  that  arguments  were  provided  to  show  that  could  be  further  reduced 
to 


°12 


Re 


iya00/a 


22  66 


* V»22'«64 


(A38) 


In  this  discussion,  it  will  be  shown  that  dropping  the  first  term  in  eqn. 
(A22)  in  order  to  obtain  eqn.  (A38)  is  essentially  equivalent  to  dropping 


9x2 


in  the  shear  strain. 
Note  that 


3u 


l 


3u. 


i/a00/a 


22  66 


(a  1 1 a 1 2 a12 
, 3/ / 3 * * 


■ 


where  eqns.  (A27)  and  (A9)  have  been  used.  Also, 


12  3x, 


Hi  e 

3x,  nr  e 


i'/a22/a66 


66  ' ^ I (*,  * vw  *2),/! 


(A45) 


where  eqn.  (A35)  has  been  used.  Note  that  eqns.  (A38)  and  (A45)  are  identical 

However,  eqn.  (A44)  cannot  be  discarded  without  a closer  examination.  The 

first  term  in  eqn.  (A44)  is  the  same  as  the  first  term  in  eqn.  (A22) , which 

has  been  dropped  in  accordance  with  foregoing  arguments.  But  the  second  term 

/ a 1 1 3 1 2 a 1 2 \ 

in  eqn.  (A44) , except  for  the  constant  | - ] , is  the  same  as 

V a66  36 V 

eqn.  (A45)  which  has  been  retained.  Thus,  in  order  to  justify  the  discarding 
of  the  second  term  in  eqn.  (A44)  also,  it  is  necessary  that 


alla12  al2 


(A46) 


From  eqns.  (Al)  and  (A4) 


rallal2 


V1 2G1 2 
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(A47) 


Thus,  the  necessary  condition  represented  by  the  inequality  (A46)  is  proved. 
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GOVERNING  EQUATIONS  FOR  MODE  I DYNAMIC  FRACTURE 


IN  ORTHOTROPIC  DCB  SPECIMENS  WITH  EXTERNAL 
DYNAMIC  INTERACT  I ONS 


ABSTRACT 


The  governing  partial  differential  equations  for  dynamic  crack 
propagation  in  double  cantilever  beam  (DCB)  unidirectional  composite 
specimens  are  derived.  Only  the  Mode  I symmetric  case  in  which  the 
loading  is  along  a principal  material  axis  and  a self-similar  crack 
propagates  along  the  in-plane  orthogonal  principal  material  axis  is 
considered.  The  natural  dynamic  interactive  boundary  conditions  between 
the  DCB  specimen  and  the  tensile  testing  machine  are  emphasized.  A 
finite-difference  solution  scheme  is  presented. 


INTRODUCTION 


In  the  present  study  of  the  dynamic  fracture  of  composite  materials, 
one  of  the  proposed  experimental  specimens  is  a double  cantilever  beam 
(DCB)  which  may  be  modeled  as  a Timoshenko  beam  on  a generalized  founda- 
tion. The  composite  material  in  the  beam  is  modeled  as  a homogeneous 
orthotropic  material  and  the  analysis  is  restricted  to  the  Mode  I symmetric 
case  in  which  the  loading  is  along  a principal  material  axis  and  a self-similar 
crack  propagates  along  the  in-plane  orthogonal  principal  material  axis. 

The  analysis  is  similar  to  that  by  Kanninen  [1]  who  derived  analogous 
governing  equations  for  a homogeneous  isotropic  DCB  subjected  to  constant 
displacement  loading. 

In  addition  to  the  governing  partial  differential  equations  of 
the  interior  of  the  beam,  the  dynamic  interactive  boundary  conditions 
between  the  DCB  test  specimen  and  the  tensile  testing  machine  are  derived. 
Kousiounelos  and  Williams  [2]  considered  the  dynamic  fracture  DCB  problem 
and  derived  the  governing  equations  but  did  not  consider  interactive 
boundary  conditions.  For  completeness,  the  derivation  of  the  governing 
partial  differential  equations  will  be  presented  here  as  well  as  the 
derivation  of  the  natural  boundary  conditions  which  account  for  the  dynamic 
interaction  between  the  test  specimen  and  the  tensile  testing  machine. 
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EQUATIONS  OF  MOTION  AND  BOUNDARY  CONDITIONS 


The  structure  to  be  analyzed  is  shown  in  Fig.  1.  From  Fig.  lb 
it  can  be  seen  that  there  are  two  regions  of  the  beam  where  different 
equations  of  motion  are  applicable.  These  are  the  cracked  region  where 
no  elastic  foundation  exists  and  the  uncracked  region  where  a generalized 
foundation  is  present.  In  deriving  the  equations  of  motion  for  both  of 
these  regions,  the  singularity  function  <x-a>°  which  is  zero  for  negative 
values  in  the  brackets  and  is  unity  otherwise  will  be  used.  A variational 
derivation  of  the  equations  of  motion  will  be  presented. 

The  variational  indicator  for  Hamilton's  principle  is  [3] 


where  t is  time,  are  the  generalized  coordinates,  and  Ej  are  the 
generalized  forces  corresponding  to  the  nonconservative  forces.  The 

•k 

kinetic  coenergy  (T  ) and  the  potential  energy  (V)  functionals  are  given 


v V (Iff  +i"ci2A(!?-*f 

+ <x-a>°  -J  (kew2  + kr^2)|dx  + ^ks  jw(0,t)  - wQtJ2 


r*+  v 


The  generalized  coordinates  are  w,  the  transverse  displacement,  and  ip, 
the  in-plane  rotation.  The  beam  has  mass  density  p,  cross-sectional  area 
A,  flexural  moment  of  inertia  I,  and  shear  coefficient  < . is  the 

extensional  modulus  in  the  1 direction  and  is  the  shear  modulus  in 
the  1-2  plane,  as  indicated  in  Fig.  1.  k£  and  kf  are  the  extensional  and 
rotational  foundation  stiffnesses,  respectively.  The  crack  tip  is  at 
x = a(t).  Because  the  crack  tip  cannot  "heal",  if  at  any  value  of  x the 
function  <x-a>°  is  zero,  this  function  cannot  become  unity  at  any  later 
time.  It  is  assumed  that  the  beam  contains  negligible  structural  damping 
and  that  no  nonconservative  forces  act  on  the  system;  thus  the  E^  are  zero. 
The  tensile  testing  machine  interactions  are  modeled  in  Fig.  2 where  an 
extensional  spring  with  spring  constant  kg  and  a translational  mass  m 
represent  the  machine's  stiffness  and  inertia,  respectively.  The  cross- 
head of  the  testing  machine  is  assumed  to  move  at  a constant  displacement 

rate,  thus  one  end  of  the  spring  k is  assumed  to  move  at  a constant 

s 

velocity  wq  . 

Substitution  of  eqns.  (2)  and  (3)  into  eqn.  (1)  and  use  of  the 
calculus  of  variations  give 


V.I 


• ■ { “!•(  “‘h  0 + ($  ■ £)  - <x'a>‘k*v] 


6w(x,t) 


f0  d*["CI  3r>  + El1  + (S-*) 


(equation  continued) 
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The  necessary  conditions  for  the  variational  indicator  in  eqn.  (4)  to 

I 


vanish  for  arbitrary  admissible  variations  of  w and  ip  are 


Eqns.  (5)  are  the  partial  differential  equations  of  motion  and  eqns.  (6) 
and  (7)  are  the  natural  boundary  conditions.  In  deriving  eqns.  (5),  (6) 


and  (7)  it  has  been  assumed  that  Gl2,  A,  I,  kg  and  kf  are  independent 
of  x and  t. 


FOUNDATION  PARAMETERS 


The  elastic  foundation  parameters  represent  the  effects  of  the 


extensional  and  bending  forces  transmitted  across  the  uncracked  plane  of 


the  specimen.  The  model  which  is  used  is  shown  in  Fig.  3.  From  Fig.  3(b) 


the  extensional  equilibrium,  compatibility  and  constitutive  relations  are 


and  E_  is  the  associated  extensional  modulus.  Combining  the  relations  in 


From  Fig.  3(c),  the  rotational  equilibrium,  compatibility  and  constitu 
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CRACK  PROPAGATION  ENERGY  CRITERION 

In  developing  the  crack  propagation  energy  criterion,  it  is 
assumed  that  heat  generation  and  heat  propagation  effects  are  negligible 
throughout  the  body  and  that  the  body  remains  linearly  elastic  except 
for  a small  undefined  crack-tip  region  in  which  the  fracture  process  is 
occurring.  So,  by  the  conservation  of  mechanical  energy,  the  incremental 
energy  absorbed  by  the  extending  crack  must  be  equal  to  the  incremental 
work  done  on  the  system  minus  the  incremental  mechanical  energy  released 
by  the  system.  Thus,  for  a given  increment  of  crack  advance  da,  this  may 
be  expressed  as 

Rb  da  = d/^-  d^ 
or 

where  R is  the  dynamic  fracture  toughness  of  the  composite,  2^is  total 
external  work  done  on  the  specimen  by  both  the  upper  and  lower  machine 
loadings,  and  $ is  the  total  mechanical  energy  in  the  system  containing 
both  halves  of  the  DCB  specimen  as  well  as  both  the  upper  and  lower  testing 
machine  elements. 

The  total  mechanical  energy  in  the  system  is 

Tt  + VT  (14) 

where  TT  and  VT  are  the  total  kinetic  and  total  potential  energies, 
respectively.  Because  the  system  is  linear 

Tt  = 2T*  and  VT  = 2V  (15) 


(12) 

(13) 
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where  T and  V are  defined  in  eqns.  (2)  and  (3).  So, 


a?  * 

d&=  9 dT_  dV 

dt  dt  1 dt  ’ 


Substituting  eqns.  (2)  and  (3)  into  eqn.  (16),  and  performing  the 
differentiation  indicated  in  eqn.  (16)  give 


(16) 


+ <x-a>°kr\|>  } dx 


+ <x-a>°k  w i dx 
e 


i 


- 2ks[w(0,t)  - wQt]  Co  + (kew2  + k^2)  I (-  |f) 


(17) 


where  the  boundary  conditions  represented  by  eqns.  (6)  and  (7)  have  been 
used  in  reducing  eqn.  (17)  to  the  form  given.  By  eqns.  (5),  the  flower 
brackets  in  eqn.  (17)  vanish,  and  thus  eqn.  (17)  reduces  to 


#-  ~2ks  - V>0  - at  (ke"’+kr*’)  x.a  • 


(18) 


The  work  done  on  the  system  is  input  via  the  testing  machine 
stiffness  by  the  constant  velocity  wq.  (Refer  to  Fig.  2.)  The  upward 
force  in  the  machine's  compliance  is  k^Jw^t  “ w (0 , t )J  . So,  the  total 


work  increment  is 


= 2k  [w  t - w(0,t)ldw 
s L o J o 


(19) 
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...  . 


■.  m.  3-.  a.  . i, 


Hr  w t* 


from  which 


2k  [w  t - w(0,t)lw 
dt  s[  o Jo 


Substitution  of  eqns.  (18)  and  (20)  into  eqn.  (13)  gives 


Rb  = ^kfiw2  + kfij/2j  j , 


Using  the  dynamic  surface  fracture  energy  y , eqn.  (21)  may  be  written 

cl 


^‘i[ke”J  + tr*2]x-a 


where  the  definition  R = 2y^  has  been  used.  Further,  these  results  can 
be  expressed  in  terms  of  the  dynamic  crack  extension  force  by  noting 


= 1 \&jT 

d ~ b Ida 


which  may  be  rewritten  as 


S'f  = £ |"k  w2  + k i^2~| 
d b [_  e r Jx=a 


where  k and  k are  given  in  eqns.  (9)  and  (11),  respectively, 
e r 


is  also  called  the  dynamic  energy  release  rate. 


w T: 


NONDIMENS IONAL  EQUATIONS 


In  preparation  for  solving  the  governing  equations  it  is  convenient 


to  nondimensionalize  the  equations  of  motion  and  the  boundary  conditions 


the  singularity  condition  <x-a>°  can  be  restated  as  <0  -9> 


Also,  it 


is  convenient  to  define  the  following  nondimensional  variables 


Substitution  of  eqns.  (26)  into  the  first  of  eqns.  (5)  gives 


Substitution  of  eqns.  (26)  into  the  second  of  eqns.  (5)  gives 


If  the  DCB  has  a constant  rectangular  cross  section  such  that 


eqns.  (27)  and  (28)  become,  respectively, 
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ax2 

and 


0 + ^(1-) 


<0  -0>° 
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5G12  y = J_  3^Y 

ei  " 12  ax* 


(31) 


Also,  if  eqns.  (26)  are  substituted  into  the  natural  boundary  conditions, 
eqns.  (6)  and  (7)  become,  respectively, 


I 


An  important  part  of  the  dynamic  fracture  computation  is  the 
evaluation  of  the  energy  and  work  expressions  during  crack  propagation. 
Substitution  of  eqns.  (26)  into  eqn.  (2)  gives 


£-^kLh[(^)2+M£)> 


3W(0,t) 

9t 


Substitution  of  eqns.  (26)  into  eqn.  (3)  gives 


Substitution  of  eqns.  (26)  into  eqn.  (19)  gives 


J 


FIN TTF.-DI  FFERENCE  SOLUTION  SCHEME 


Equations  (30)  and  (31)  may  be  written  in  their  respective 
finite-difference  formulations  as 


— — |w(C+AC,  t)  - 2W(C,T)  + W(C-AC,  t)1 

(AO2  ' 


12e 

- 2^  {y(OA?,  t)  - Y(C-AC,  T)|  - <ec-0>°  5^  W(;,t) 
= - |w(C,  t+At)  - 2W(C,t)  + W(C,t-At) j 

i \ ^ ' 


10G12(At)‘ 


— {y(c+AC,  t)  - 2Y(c,x)  + Y(c-A C,  t)} 
(AC)2  1 1 


|w(C+AC,  T)  - W(C-AC,  T)|  - [2  + <0c-0>oJy(C,' 


{ Y(C,  T+At)  - 2Y(C,T)  + Y (C,  t-At ) | 

12 (At) 2 1 ' 


where  A is  the  incremental  finite-difference  operator.  Solving  eqn.  (37) 
for  W(c,  t+At)  gives 


w(c,  t+At)  = 


kfal  1 i 


10G 


E^(f^)  {W(^,  T)  + W(C-AC.  T)} 


E1 

5G12 

(At)2 

E1 

AC 

|20G12 

/AiY 

X^Cj 


(Ai)2<ec-e>°  w(?lT) 


- W(C,  t-At), 


(39) 


Solving  eqn.  (38)  for  Y(c,  t+At)  gives 


Y(c,  t+At)  = 12  (|^)2{y(C+AC,  t)  + Y(C-AC,  t)} 


, 60G12  (At)2 


Ex  A? 


|w(C+AC,  t)  - W(C-AC,  t)| 


{24(ac  ) - 2 + 60  ^ (Ax>2[2  + <ec“e>° J|y (C ,t) 


- Y(C,  t-At). 
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Eqns.  (39)  and  (40)  are  the  recursion  relations  to  be  used  in  the  computa- 
tions. The  associated  boundary  conditions  may  be  written  by  expressing 
eqns.  (32)  and  (33)  in  their  respective  finite-difference  formulations  as 


Y(-AC,  T)  = Y(AC,  T) 


W(0,t+At)  = g2  \!kG12A  h ~ ,T>  - Y (0 , X )1 

m 12  hL  2AC  j 


2AC 


- kg[w(0,T)  - Wqt]  - W(0,t-At) 
+ 2W(0,t) 


at  C=0  (41) 


and 


H(s  + 4«-T)-“(h'4c’T)  + 2A?Y(h’T) 


“ ' "h 


Further,  the  initial  conditions  can  be  obtained  by  solving  the  static 
problem  of  a Timoshenko  beam  on  an  elastic  foundation  for  the  loading 
conditions  at  x=0. 
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(42) 
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J 


From  eqn.  (34)  and  by  using  backward  differences  for  the  time 


derivatives,  the  kinetic  coenergy  may  be  written  as 


2T 

RA 


1 E1 


24  E 


E 


2 Over  Beam 
Length 


3W(C.t)  - 4W(S,  t-At)  + w(c.  T-2AT) 


2At 


12 


3Y(C,t)  - 4Y(C.  t-At)  + Y(£,  t-2At ) 
2At 


AC 


_1 1 m 

12  E2  phA 


3W(0,t)  - 4W(0.  t-At)  + W(0,  t-2At) 


2At 


(43) 


Similarly,  from  eqn.  (35)  and  by  using  central  differences  for  the  spatial 
derivatives,  the  potential  energy  may  be  written  as 


2V 

RA 


E 


Jl  fl 

_ „ . 24  E, 

Over  Beam  I 2 

Length 


Y ( C+AC . T ) - Y(C~AC,  t) 
2AC 


_5_  fl2 
12  E„ 


+ <0  -0> 
c 


W(C+AC,  T)  - W(C-AC,  T)  _ y(c>t) 
2AC 


[w2(c. 


l»x)j 


k 


2E2b 


[w(0,t)-Wqt]  . 


(44) 


The  corresponding  finite-difference  equation  for  the  work  expression  is 
identical  to  eqn.  (36).  Further,  it  is  noted  that  the  forms  of  eqns.  (43) 


and  (44)  are  written  differently  so  as  to  ensure  that  the  approximate 
errors  in  each  are  of  the  same  order. 


It  is  necessary  to  establish  an  integration  step  size  in  order  to 
perform  the  above  finite-difference  calculations.  A thorough  error 
analysis  of  the  governing  differential  equations  [that  is,  eqns.  (30) 
and  (31)]  provides  the  following  restrictions  [4]: 


AC  - \/  10G- 


< /-L 

Ac  -\/l2 


where  the  appropriate  ^ is  the  smaller  of  the  two. 
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Models  for  calculating  the  foundation  parameters  of  the  orthotropic 
double  cantilever  beam. 


SECTION  h 


NUMERICAL  STABILITY  AND  CONVERGENCE  CRITERIA 
FOR  EQUATIONS  FOR  DYNAMIC  FRACTURE  IN 
DCB  COMPOSITE  SPECIMENS 


ABSTRACT 

The  stability  and  convergence  criteria  for  numerical  solutions 
of  the  governing  equations  for  dynamic  crack  propagation  in  DCB  uni- 
directional composite  specimens  are  determined.  The  stability  criteria 
are  given  explicitly  and  it  is  concluded  that  the  convergence  require- 
ments will  be  met  automat ical ly  by  solutions  which  satisfy  the  stability 


■ — I .ill ..HI  M1..I ■ mm 


INTRODUCTION 


The  governing  differential  equations  for  dynamic  crack  propagation 
in  double  cantilever  beam  (DCB)  specimens  generally  require  a numerical 
solution.  These  equations  are  frequently  restated  in  their  finite- 
difference  approximations  as  in  the  cases  for  the  homogeneous  isotropic 
DCB  [1]  and  for  the  homogeneous  orthotropic  DCB  [2].  The  numerical  evalu- 
ation of  the  finite-difference  equations  requires  the  selection  of  temporal 
and  spatial  interval  sizes. 

Heuristic  suitable  interval  sizes  for  the  finite-difference  isotropic 
DCB  equations  have  been  obtained  by  trial-and-error  [1].  However,  no 
analytical  verification  of  suitable  interval  sizes  has  been  obtained  for 
this  problem.  Furthermore,  suitable  interval  sizes  for  the  orthotropic 
DCB  have  not  been  obtained  either  by  trial-and-error  or  analytically.  The 
purpose  of  this  report  is  to  present  an  analytical  derivation  of  suitable 
interval  sizes  for  the  solution  of  the  finite-difference  orthotropic  DCB 
equations  for  Mode  I dynamic  crack  propagation.  The  analogous  suitable 
interval  sizes  for  the  isotropic  DCB  are  contained  within  these  results. 

If  5^  represent  the  exact  solutions  of  the  partial  differential 

equations  and  5j  represent  the  exact  solutions  of  the  finite-difference 

equations  and  5^  represent  the  numerical  solutions  of  the  finite-difference 

equations,  then  determining  the  conditions  under  which  (H)  -+•  5^1,  , , „ 

\ j i/i«j*l,2, . . . ,N 

is  the  problem  of  convergence  and  determining  the  conditions  under  which 
IH’  -*■  5")  , . is  small  throughout  the  region  of  integration  is 

the  problem  of  stability,  where  N is  the  number  of  differential  equations. 

In  the  analyses  which  follow,  the  problem  of  stability  will  be  emphasized 


I 
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as  it  will  be  shown  that  for  the  honogeneous  orthotropie  DCB  equations 
solution  stability  will  ensure  solution  convergence. 


THEORETICAL  ANALYSIS 


Governing  Equations 


The  governing  differential  equations  for  Mode  I dynamic  crack 
propagation  in  orthotropic  DCB  specimens  are  [2] 


i 


32y  +lzi2 

^2  E, 


(fM- 


<0  -6>°  — =—  Y = J2  — 
C E1  12  3t2 


sfw  3Y  fl>0  ^2 

V6  5G12 


J12  3t2 


These  equations  are  written  in  terms  of  the  following  nondimens ional 


variables: 


Independent  - 


Dependent 


- "‘“ft)' 


Y-  *h(^) 


where 


= longitudinal  coordinate  along  the  beam  (crack 
propagation  direction) 


transverse  displacement  of  beam  on  elastic  foundation 
rotation  of  beam  on  elastic  foundation 
half-height  of  DCB 


! 


' I 


I 


} 

I 


t 


width  of  DCB 

mass  density  of  the  DCB 

extensional  modulus  in  x-direction 

extensional  modulus  transverse  to  x-direction 

1-2  plane  shear  modulus 
2E2b 

h 

KG12bh 

2~ 


shear  coefficient 

K w2  + K lii2 
e r 

2b"d 

dynamic  surface  fracture  energy 

singularity  function  which  is  one  (1)  for  positive 
argument  and  zero  (0)  otherwise. 


The  respective  finite-difference  formulations  of  eqns.  (1)  and  (2)  are 


~ {y(C+A 5,  t)  - 2Y(?,t)  + Y(C-AC,  x)J 
(AO2  ' > 

+ E^ACT  {w^+a?*  ~ W(C-AC,  T)}  - ~ [2  + <0c-6>°]y(c,t) 

" {Y(5.  t+At)  - 2Y(C,t)  + Y(c,  t-At)1  . (3) 

12 (Ax) 2 1 f 
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--hw- 


— — fw(C+AC,  T)  - 2W(?,t)  + W(£-AC,  T)i 

(AO*  1 ’ 

12e 

- 2fe{Y<^.  T>  - Y<^.  t)|  - <ec-e>° 


|w(ct  t+At)  - 2W(5,T)  + W(C,  T-AT)i 

a-r\2  ‘ ' 


iog12(At) 


where  A is  the  incremental  finite-difference  operator. 


Error  Equations 


It  is  useful  to  define  the  following  variables: 


yj>kSY(c,T)  ; Vttl  k±1  = V«±4C,  t±4t) 


UJ,k  = = MJ±l.k±l  5 Ti4T> 


where  j and  k are  integers.  If  Y , and  W.  , represent  the  exact 

Ji11  3 >K 

solution  of  the  finite-difference  equations  (3)  and  (4),  then  the 
numerical  solution  of  the  finite-difference  equations  may  be  written  as 


(Y+A)j  k and 


where  A and  y are  numerical  errors  associated  with  Y and  W,  respectively. 
Substitution  of  eqns.  (5)  into  eqns.  (3)  and  (4)  gives 
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Substracting  eqn.  (7)  from  eqn.  (9)  and  eqn.  (8)  from  eqn.  (10)  gives 


Vi.k  ~ 's.k 

<«): 2 4C 


- (A+B ) X 


j,k 


„ ~ 2X1,k  + X1,kH 

(At)2 


(ID 


and 


a V-iA : 2uj,k i . A . c 

(AC)2  AC  j,k 


Uj,k+1  ~ 2yj,k  + V,i,k-1 

(At)2 


(12) 


It  can  be  seen  that  the  errors  X and  p must  satisfy  the  same  equations  as 
their  respective  associated  nondimensional  variables  Y and  W.  This  is  a 
consequence  of  the  linearity  of  the  governing  partial  differential  equations 
(1)  and  (2). 


Solution  Stability 

The  solution  stability  analysis  which  fo'  .ows  is  generally  attri- 
buted to  J.  von  Neumann  although  the  first  detai  ed  discussion  was  pre- 
sented by  O'Brien,  Hyman  and  Kaplan  [3].  The  procedure  is  to  form  Fourier 
expansions  of  lines  of  errors  and  to  follow  the  progress  of  the  general 
terms  of  the  expansions.  For  the  problem  under  consideration,  these  Fourier 
expansions  may  be  written  as 
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.... 


■■  J*. 


e e 
n 


(13) 


r' 


i r 


c.  i 


_ 18  5 18  5 

1(5)  = ZJ  6 e n and  p(5)  = £)  e e ° 


At  each  line  of  computation, groups  of  errors  are  introduced  which  propagate 
forward  in  accordance  with  eqns.  (11)  and  (12). 


Because  of  superposition,  it  is  necessary  to  consider  only  the 


185 


single  error  term  e where  8 is  any  member  of  the  |$n| , each  of  which  is 

185 


any  real  number.  So,  solutions  of  eqns.  (11)  and  (12)  which  reduce  to  e 
when  T=0  are  sought.  In  accordance  with  [3],  let 


e“Tei65 


(14) 


where  a = a (8)  is  generally  complex.  In  order  that  the  original  error 
185 


shall  not  grow  as  x increases,  it  is  necessary  and  sufficient  that 


1,131, 


Substitution  of  eqn.  (14)  into  eqn.  (11)  and  eqn.  (14)  into  eqn.  (12) 
gives, respectively. 


2(cos8A5-l)r2-  D(n-2+  ^ ) - (A+B)Ax2  A(eigAC  - l)rAx 


-A(ei6A^  - l)rAx 


2A(cos8A5~l)r‘ 


- (n-2+  " CAx2 


(15) 


- Al 
A5 


where  n = e ‘ and  r = -77  . A non-trivial  solution  requires  that  the 


_ aAx 
= e 

determinant  of  the  coefficient  matrix  be  zero.  Thus, 
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[*-ll'l  + [A2(ei6A?-l)r2  - (A+B)'!'  - C4>]Ax2  + [(A+B)C]Ax4  *=  0 (16) 


where  ♦ • 2(cos0AC-l)r2  - D(n-2+  — ) 


and  V ■ 2A(cos0A£  - l)r2  - (n-2+  — ) 

n 


Because  the  coefficients  of  Axn  are  of  the  same  order  and  by  invoking  the 
requirement  that  Ax  « 1,  it  can  be  seen  that  the  first-order  approximation 
to  the  solution  of  eqn.  (16)  is  satisfied  by  requiring  that  the  coefficient 
of  Ax°  is  approximately  zero.  Thus,  $>•4'  = 0 gives 


A2  - 2(cos0Ac-l)(A  + i)  r2A  + -y  (cos0A£  - l)2r4  = 0 (17) 


where  A = n-2+—  • The  solution  of  eqn.  (17)  may  be  obtained  by  use  of 
the  quadratic  equation  as 


I 

I 

I 

I 

I 

I 

I 


A = (cosgAC  - l)r2 


(a^)±V(a  + 5)2  -T 


(18) 


Noting  that  (a+q)  “ ^ = (A_d)  ’ eqn*  gives 


IN 

|lj  2(cos0AC  -1  )r2  . 


(19) 


10G 


Note  that  because  A = 


12  1 

— and  D = are  always  real  and  positive 


and  because  -1  cosBAC  1,  the  maximum  value  of  both  A^  and  is 
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.....  ■ . i .. 


; dt,  ♦ , 


; < 


zero.  Also,  the  minimum  values  of  A and  A2  occur  when  cosgA^  = -1  . 
From  the  definition  of  A, 


ai,2  ‘ "1,2  - 2 + n 


1,2 


or. 


ni,2  ~ (2  + Al,2)ni,2  + 1 = 0 » 


(20) 


which  has  the  solutions 


*1,2  = Yl,2  - 1 


(21) 


1 2 

where  Yj.  2 = 1 + 


^ • Recall  that  the  necessary  and  sufficient  con- 
oAt| 


I Pi  2 £ 1 • It  can  be  shown,  that  these 


ditions  on  n are  that  |e' 
necessary  and  sufficient  conditions  are  satisfied  if  ^ 2|  < 1.  Note  that 

lf  |Y1,2|  - 1 ’ ^yi,2  ~ 1 = = i/l-Y^2  where  1 H VCT-  Thus, 

from  eqn.  (21)  it  can  be  seen  that  2 are  complex  numbers  with  a modulus 

of  one  (1).  Further,  the  requirements  on  y may  be  written  as 

1 , z 

A, 


1 2 

-1  < 1 + < 1 


(22) 


where  the  definition  of  Y^  2 has  been  used.  Recall  that  the  maximum  value 

°f  Al,2  is  zero;  thus,  the  right-hand  inequality  has  been  satisfied  already, 

The  left-hand  inequality  imposes  the  other  limit  on  A , namely 

1,2 


Al,2  i -4  • 


(23) 


J 
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Finally,  by  noting  that  the  minimum  values  of  A are  obtained 
when  cosBAC  = -1,  eqns.  (19)  and  (23)  may  be  combined  to  give  the  simul- 
taneous requirements 


2A  (-2)  (r2)  > -4 


2 ^ (-2) (r2)  > -4 


r < D 


By  inserting  the  definitions  for  r,A  and  D into  eqns.  (24),  the  integration 
interval  size  requirements  become 


< 

AC  - 


¥-  < 

AC  - 


where  the  appropriate  is  the  smaller  of  the  two.  Thus,  eqns.  (25) 

represent  the  criteria  to  ensure  solution  stability,  Furthermore,  eqns. 
(25)  require  that  independent  of  the  material  constitutive  parameters 
and  G ^ » r should  always  be  less  than  ^1/12  . Therefore,  eqns.  (25) 
contain  within  them  the  solution  convergence  criteria  also  because  r<l 
is  the  necessary  and  sufficient .condition  for  solution  convergence  [3,4]. 
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SECTION  5 

CONDUCTIVE  GRID  VELOCITY  MEASUREMENT  SYSTEM 


ABSTRACT 


A conductive  grid  system  which  is  capable  of  measuring  the 
crack-tip  velocity  profile  has  been  designed,  built  and  proof-tested 


VELOCITY  MEASUREMENT  SYSTEM 


A conductive  grid  velocity  measurement  system  has  been  designed, 
built  and  proof-tested.  This  system  which  Is  capable  of  providing  the 
crack-tip  velocity  profile  (crack-tip  velocity  as  a function  of  crack 
length)  is  shown  in  Fig.  1.  Lines  of  brittle  conductive  silver  paint 
are  painted  perpendicular  to  the  anticipated  crack  path.  The  silver 
paint  grid  forms  part  of  an  electrical  circuit  whose  elements  are 
sequentially  interrupted  as  the  crack  advances  and  fractures  the  brittle 
paint  lines.  The  crack  velocity  can  be  calculated  from  knowledge  of  the 
location  of  the  conductive  lines  and  the  time  at  which  each  line  fractures. 

The  design  of  the  circuit  and  the  choice  of  the  resistors  are 
such  that  the  output  DC  voltage  decreases  by  an  equal  amount  each  time 
the  fracture  of  a silver  paint  conductive  line  occurs.  Because  the 
graphite  epoxy  specimens  are  conductive,  a thin  insulating  coating  of 
epoxy  is  applied  to  the  specimen  surface  prior  to  the  application  of  the 
silver  paint  grid.  Preliminary  tests  have  shown  that  the  fracture  strains 
of  the  silver  paint  and  the  epoxy  matrix  are  essentially  the  same. 

A modified  video  tape  recorder  is  used  to  record  the  output  of 
the  circuit.  The  use  of  the  tape  recorder  eliminates  any  trigger 
requirements  in  capturing  the  dynamic  circuit  output  during  fracture 
propagation.  The  recorder  can  be  started  well  in  advance  of  the  rapid 
fracture  event  and  run  throughout  the  loading  of  the  specimen.  The 
video  tape  recorder  has  a gated  playback  capability  which  allows  the 
selection  of  any  portion  of  the  tape  to  be  played  in  the  "still"  mode. 

The  tape  recorder  output  is  displayed  in  the  "still"  mode  on  an 


Unidirectional  90°  (fibers  perpendicular  to  the  loading  axis) 
Hercules  AS/3501  graphite  epoxy  panels,  35-56  cm  wide  x 1 5 . 24  cm  high 
x 0.127  cm  thick  (14  in  x 6 in  x 0.05  in)  with  a 5-08  cm  (2  in)  Mode  I 
edge  notch  have  been  tested.  By  varying  the  loading  conditions  in  the 
tensile  test,  crack  velocities  from  30.48  m/sec  (100  ft/sec)  to  2743  m/sec 
(9000  ft/sec)  have  been  measured. 


I 


Experimental  circuit  for  crack  velocity  measurement  using  the  conductive  grid. 


SECTION  6 


CONCLUSIONS  AND  RECOHHENDAT I ONS 

The  results  obtained  so  far  have  been  described  in  Sections  1 
through  5 of  this  report  and  are  summarized  in  the  abstracts  of  the 
respective  sections.  Therefore,  these  results  will  not  be  repeated 
here  except  for  the  brief  comments  in  the  next  paragraph. 

In  formulating  a total  framework  for  dynamic  fracture,  three 
major  categories  of  the  problem  may  be  defined.  These  are 

1.  Intralaminar- interfiber  fracture, 

2.  Intralaminar-transfiber  fracture,  and 

3.  Interlaminar  fracture. 

Note  that  Section  2 dealt  with  the  theoretical  dynamic  crack  extension 
force  for  intralaminar-transfiber  fracture.  Sections  3 and  1*  dealt  with 
the  analytical  aspects  of  both  the  intralaminar-interfiber  and  the  intra- 
laminar-transfiber fracture  problems.  Section  5 described  an  experi- 
mental system  which  is  capable  of  measuring  crack-tip  velocity  profiles 
in  the  determination  of  experimental  dynamic  fracture  toughness. 

The  ultimate  goals  of  this  research  are  (1)  to  establish  firm 
theoretical  and  experimental  bases  for  dynamic  fracture  and  crack  arrest 
in  fiber  composites,  and  (2)  to  devise  standardization  techniques  and 
procedures  for  the  measurement  of  dynamic  toughness  and  crack  arrest 
capability  of  advanced  fiber  composites.  At  this  point,  there  remains 
a considerable  number  of  specific  aspects  of  dynamic  fracture  and  crack 
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arrest  in  composites  which  require  both  experimental  and  theoretical 
analysis.  A detailed  delineation  of  the  various  research  requirements 
to  achieve  these  goals  is  beyond  the  scope  of  this  program.  Therefore, 
a limited  list  of  specific  recommendations  are  offered: 

1.  Experimentally  measure  the  dynamic  fracture  toughness 

of  Hercules  AS/3501  graphite  epoxy  composites  for  intra- 
laminar-interfiber fracture. 

2.  Conduct  numerical  analysis  for  the  orthotropic  DCB  indi- 
cated in  Sections  3 and  A of  this  report. 

3.  Derive  a dynamic  crack  extension  force  via  the  velocity 
correction  factor  procedure  for  intralaminar-interfiber 
fracture. 

A.  Establish  the  required  correlations  to  bring  the  three  (3) 
previous  recommendations  into  consistency. 

These  four  (A)  recommendations  provide  a meaningful  framework  for  the 
intralaminar-interfiber  dynamic  fracture  problem  which  based  upon  this 
report  can  be  conducted  in  a straightforward  manner. 
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